Constraints from T c and the isotope effect for MgB 2 
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With the constraint that T c = 39 K, as observed for MgB2, we use the Eliashberg equations 
to compute possible allowed values of the isotope coefficient, /3. We find that while the observed 
value /3 bs = 0.32 ± O.Of can be obtained in principle, it is difficult to reconcile a recently calculated 
spectral function with such a low observed value. 
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I. INTRODUCTION 

The discovery of an isotope effect in superconducting 
Hg (T]J^] had a significant influence on the 

subsequent development of the theory of superconduc- 
tivity. Since that time the isotope effect has been mea- 
sured for almost every superconductor where possible, 
with results that span the range of complete agreement 
with BCS expectations to a null result ||. In this latter 
instance, the implications for mechanism have differed, 
depending 

on the circumstance. For example if T c is low, argu- 
ments based on competition of the electron-phonon in- 
duced attractive coupling with the repulsive Coulomb 
pseudopotential suggest that the isotope coefficient 
should be close to zero in this case . 

Qualitatively, the argument is as follows: increasing 
the phonon frequency (by reducing the ionic mass) serves 
to increase T c because the energy scale that governs T c 
is the phonon energy. On the other hand T c is also 
largely determined by the ability of phonons to medi- 
ate a retarded attractive interaction that avoids the in- 
stantaneous Coulomb repulsive potential. This requires 
a large difference between the phonon energy scale (ap- 
proximately the Debye frequency, uir>) and the energy 
scale of the Coulomb potential (roughly half the band- 
width, D/2). Raising the phonon energy scale (through 
a lower mass isotopic substitution) brings these two en- 
ergy scales closer together, and therefore acts to reduce 
T c . This effect offsets the gain in T c expected due to the 
first argument, and can lead, in some cases to an isotope 
coefficient close to zero. 

The discovery of superconductivity in MgB2 at 39 K 
f| has naturally prompted questions about mechanism, 
fsotopic substitution of both Mg and B has resulted in 
an overall isotope coefficient of 0.32 ± 0.01 As- 
suming the principle of superposition holds for the iso- 
topic substitutions, in the absence of Coulomb repulsion, 
a value of (3 = 0.5 is expected. Here, j3 = where 
f3i = —dh\T c /dh\Mi, and Mj is the mass of the i th ion, 
and the summation is over all the ions in the compound 
H . The simplest version of Eliashberg theory |l(J then 
requires a Coulomb pseudopotential to explain this re- 



duced value. 

In fact, as explained there, and more recently in Ref. 
|^ , one can gain a qualitative understanding of this re- 
duction through the modified McMillan equation pLQ , 
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Here A is the electron-phonon coupling constant, u>\ n is 
the logarithmically weighted average phonon frequency, 
and /x* is the Coulomb pseudopotential. It is important 
to note, for comparisons which we will make later, that 
the McMillan equation was derived with fj,* = /x*(cJi n ). 
That is, the cutoff associated with the pseudopotential is 
the characteristic phonon frequency itself. 

The approximate formula for obtaining the pseudopo- 
tential at any frequency Wj, given the 'bare' potential fi 
associated with the bandwidth energy scale, D/2, is 



fi(D/2) 



l + /z(-D/2)ln 



D/2 



(2) 



Equation ^| requires assumptions of particle-hole symme- 
try and the smallness of the parameter T c j Ui . Removal of 
these assumptions is further discussed in Ref. Jll|. Most 
importantly, however, for this relation to work requires 
Wi^6u;i n . Nonetheless, for purposes of comparing nu- 
merical solutions of Eliashberg theory to the McMillan 



equation, we will apply Eq. 2 

Based on Eqs. 
given by 



for uji — 



1|) and (|D the isotope coefficient is 



1 



1.04(1 + A)(1 



0.62A) a 



[A-m*(1 + 0.62A)] 
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These equations have been used previously to constrain 
A and fi* |l0|]. Very recently Hinks et al. || have ar- 
gued that their observed isotope effect is incompatible 
with a small value of //*, as sug gest ed by calculations of 
the electron-phonon parameters |1J] . The purpose of this 
paper is to examine these conclusions with numerical so- 
lutions to the Eliashberg equations. In this way not only 
can one examine the accuracy of the McMillan equation, 
but also the dependency on details of the electron-phonon 
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spectral function, a 2 F(uj) (recall that the fitted parame- 
ters in the McMillan equation were initially based on the 
Nb spectrum). In the following section we will first exam- 
ine solutions for an Einstein function spectrum, for which 
parameters can be varied simply and systematically. We 
find, upon comparison with the results from the McMil- 
lan equation, reasonably good agreement. In Section III 
we examine first the effect of broadening this spectrum, 
and then the effect of using a more realistic spectrum, 
namely the one which has been calculated with ab initio 
methods p"5| ] . A summary is provided in the final section. 

II. EINSTEIN SPECTRUM 

Inspection of the calculated spectrum in Ref. jl5| shows 
a spectrum that is dominated by one peak at approx- 
imately 75 meV. While we will later explore the con- 
sequences of this spectrum, for now we simply model 
a 2 F(uj) with an Einstein spectrum: 

a 2 F{uj) = ^5{lu - fl). (4) 

Based on the known values of T c (39 K) and (3 (0.32), 
then, for a given phonon frequency and value of [i, one 
can determine the electron-phonon coupling strength, 
characterized by A. 

Because we will vary the phonon spectrum character- 
istics in what follows (in particular the frequency scale), 
we choose to characterize the direct Coulomb potential by 
its value at some electronic energy scale associated with 
the bandwidth we choose this energy scale to be 
D/2 = 5 cV, to correspond roughly with the bandwidth 
of MgB2 fjj],[l8| • In Fig. 1 we show the results of this ex- 
ercise. In Fig. la we first show vs. fi(D/2 = 5eV), 
for a variety of Einstein frequencies. These curves are 
calculated directly from Eq. (0); note, however, that 

a different cutoff frequency has been used for each 
curve, corresponding to the phonon frequency for that 
particular case. It means that as fx(D/2) increases (in 
principle without bound), /i*(f2)saturates, and in fact, 
can barely exceed a value 0.2 [Ell. Note that for /i ss 2, 
then, with D « 10 eV, the effective Hubbard U is already 
approximately 20 eV. 

In Fig. lb the values of A required to yield T c = 39 
K are shown as a function of /i(D/2). Note that as 
/i (henceforth we will denote fi(D/2) by simply fi) in- 
creases A also saturates, specifically because, as shown in 
Fig. la, /i*(n) saturates. Moreover, the saturated values 
vary with frequency, but for the most part are reasonably 
moderate, of order unity ]l9|| . 

Finally, in Fig. lc we show the calculated isotope coef- 
ficient, /3, as a function of fi. These of course begin at 0.5 
(for fi = 0), and fall gradually, while also showing some 
signs of saturation (again, as a consequence of Fig. la). 
This figure restricts the possible parameters; we obtain 
an isotope coefficient close to the observed result with a 
sufficiently high phonon frequency. 



Thus, the observed isotope shift in MgB2 can readily 
be obtained within the framework of Eliashberg theory. 
Visual inspection of Fig. 1 shows that the value of 
must be 'high', i.e. it must exceed the 'canonical' value 
of 0.1 (as is often adopted in studies of T c ). On the other 
hand, the value required is not excessive, i.e. /Li* (O) « 0.2 
is sufficient. Moreover, the required value of A is mod- 
erate, as already commented on. The last restriction is 
the phonon frequency itself. Here, we require fairly high 
phonon frequencies, with » 100 meV. This latter con- 
straint poses two difficulties. First, the required value 
must exceed that calculated |H| by a small amount (in 
addition the coupling strength here exceeds that calcu- 
lated by a significant amount as well). Secondly, even if 
we disregard the ab initio calculation the possibility of 
high frequency phonons coupling to electrons with Awl 
has been questioned in the past H). We are unable to 
resolve the issue of the last comment; however, before we 
suggest that a discrepancy exist with the ab initio work, 
we should further investigate the impact of spectral func- 
tion shape on our results. We address this in the next 
section. 

Before we leave this section, however, it is of interest to 
compare our results with those obtained by the McMillan 
equations (|l|) and (||), as used in Ref. [||]. In Fig. 2 we 
plot (a) A, and (b) f3 vs. /i, for a few select Einstein 
frequencies. For the McMillan equations, we have used 
Eq. (J2J) (with LOi = fi) to determine /x*(f2). It is clear 
that the McMillan equation is excellent in this regime; 
the value of A is accurate to within 5%, while the value of 
(3 is somewhat over-estimated by the McMillan equation, 
but follows the trends given by the numerical solutions 
to the Eliashberg equations fairly closely. 

Nonetheless, because it is important to explore spec- 
tral shape dependence, we continue with our use of the 
numerical Eliashberg equations |hJ . 

III. LORENTZIAN SPECTRUM, AND AB INITIO 
SPECTRUM 

We begin first with some systematics. To understand 
how a realistic spectrum impacts our conclusions, we 
adopt the following straightforward procedure. First, we 
continuously broaden the Einstein mode at a given fre- 
quency into a Lorentzian with mean frequency O and 
half- width a (see Fig. 3a). We begin at a particular 
point on the curves in Fig. 1. For example, if we take 
fx = 1.2 and Q = 75 meV, then A = 1.08 and f3 = 0.355. 
Now we imagine that /x is fixed, and upon broadening, 
the value of A is readjusted (in this case it is increased) 
to maintain T c = 39 K. We now calculate (3 (as before, 
except now we are required to shift an entire spectrum) 
to see if a broadening of the spectral function decreases 
or increases (3. 

The solid curve in Fig. 3c shows our results for the pa- 
rameter regime used above (/1 = 1.2 and fi = 75 meV). 
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For calculations the values of sigma were chosen to be 
0.1, 2, 6, and 40. Clearly the effect of spectral broad- 
ening is small, and, in this regime, works to raise the 
isotope coefficient. We thus anticipate that adopting a 
more realistic spectrum will take us further from the ex- 
perimental value. Since the mean frequency is rather 
high, we also broaden the spectrum asymmetrically, as 
indicated in Fig. 3b. The result is that (3 is increased 
even further, as shown by the dashed line in Fig. 3c. 
While we have no proof, we suspect that for a given T c 
and /x, an Einstein spectrum represents a lower bound 
for the value of (3. This allows us to use Fig. 1 to make 
definitive statements concerning the possible value of /3, 
obtained for a given T c and /i (or, alternatively, A). 

We take as a most plausible spectrum for MgB2 the one 
calculated in Ref. jl5|] . Using the absolute value of the 
coupling A = 0.836 calculated from this spectrum will 
lead to an isotope coefficient very close to 0.5, in clear 
disagreement with experiment [|22f . In fact this can be 
easily estimated on the basis of Fig. 1. With a value of 
fi*(uj c = leV) ~ 0.1, then /i « 0.12, and so Fig. 1 gives 
/3^0.46, regardless of what the characteristic frequency 
of the spectrum is. 

To reconcile this discrepancy, we adjust the overall 
coupling strength (maintaining their calculated spectral 
shape) in an attempt to achieve (3 ~ 0.32 with minimal 
change. As Fig. 1 shows, we will need to increase /i, and 
hence A. Thus, we assume that the calculated spectral 
strength has somehow been underestimated. The calcu- 
lated dependence of j3 on /x is shown by the dashed line 
in Fig. 4. We find that (3 remains always above 0.355 for 
H < 2 (i.e. C/<20 eV), and, for more moderate values 
of /x, is greater than 0.4. 



IV. SUMMARY 

Fig. 1 illustrates that in principle, an isotope coeffi- 
cient (3 « 0.32 with T c = 39 K is possible with a modest 
value of A. However, this does not take into account the 
objection raised long ago (20) to the possibility of a sys- 
tem to maintain high phonon frequencies and moderately 
strong electron-phonon coupling. The use of a more re- 
alistic spectrum reinforces this difficulty. The value of (3 
is raised, and so a spectrum with an even higher char- 
acteristic phonon frequency is required. Finally, using 
the spectral shape calculated in Ref. jl5| yields an even 
higher value of /3, so that a value of (3 « 0.4 is obtained, in 
clear disagreement with experiment. Possible solutions to 
this discrepancy might include the effect of the band edge 
about 0.5 eV from the Fermi energy. As discussed in the 
second entry of Ref. Jr| a low hole concentration leads to 
a reduced isotope effect. However, for MgB 2 the energy 
difference between the band edge and the Fermi level 
is still about 5-10 times the characteristic phonon fre- 
quency. A second possibility is the role of anharmonicity 
p3[ in the primary phonon which couples to the electron. 



Unfortunately, no formalism presently exists that allows 
us to compute T c from an electron-phonon interaction 
in the case where the phonons are anharmonic. Finally, 
there remains the possibility that the mechanism is not 
the electron-phonon interaction p^ j. Some enlightening 
data would be a measurement of the isotope coefficient 
for either hole or electron doped MgB2. The primary ef- 
fect would be T c (n) itself, as has been discussed in Refs. 
p^ , p5[ , but a measurement of f3(n) would also help differ- 
entiate between doping effects due to the electron-phonon 
interaction, and doping effects due to some other mech- 
anism. 
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FIG. 1. a) The pseudopotential fi*(u c = Q), b) electron-phonon coupling constant A, and c) isotope coefficient /3, as a 
function of [i. The value of fi is varied while T c is held fixed at 39 K. These curves represent solutions to the linearized 
Eliashberg equations with a cutoff lo c = 5 eV, with an Einstein spectrum with a frequency as shown. Smaller values of f3 are 
obtained only for spectra with high frequencies, and, in that case, with high values of [i. 
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FIG. 2. Comparison of the value of a) A, and b) /3, obtained with the McMillan equation compared to a full Eliashberg 
solution for an Einstein spectrum, as in Fig. 1. For all values of fx the agreement is excellent for A, and qualitatively very good 
for 13. 
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FIG. 3. Dependence of isotope coefficient (3 on the broadening parameter a. In a) and b) the actual spectra used are drawn 
(see text for detail). In c) the value of (3 is plotted as a function of a. Clearly, (3 increases with broadening (while T c and n are 
held fixed). 
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FIG. 4. (3 as a function of /i, as determined by scaling the calculated spectrum of Kong et dl. ([LB) to keep T c — 39 K. To 
achieve f3 < 0.4 requires a very large value of fi. The experimentally observed value of /3 = 0.32 ± 0.01 is very difficult to obtain 
with this calculated spectrum. Results for the Einstein spectrum with fi = 60 and 70 meV are also shown for comparison. 
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